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j^— N ' The Schwinger-Dyson equation for the quark in the rainbow approximation at finite temperature 

("■^ I (T) is solved numerically without introducing any ansatz for the dressed quark propagator. The 

rvQ ■ dymanical quark mass-function and the wave-function renormalization are found to have non-trivial 

dependence on three-momentum, Matsubara-frequency and temperature. The critical temperature 
^ , of the chiral phase transition {Tc) and the T-dependence of the quark condensate are highly affected 

by the wave-function renormalization. We found that Tc — 155 MeV which is consistent with the 
C^ ' result of the finite temperature lattice QCD simulation. It is also found that the system is not a gas 

of free quarks but a highly interacting system of quarks and gluons even in the chirally symmetric 

phase. 
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I. INTRODUCTION 



t~^ I The vacuum of the quantum chromodynamics (QCD) is believed to undergo a phase transition to the chirally 

^D ' symmetric and deconfinement phase at high temperature and/or high baryon-density due to the asymptotic freedom 

and the plasma screening. Such a new state of matter is called the quark-gluon plasma (QGP) and is expected to be 

^^ . produced in the on-going Relativistic Heavy-Ion Collider (RHIC) at BNL and in the future Large Hadron Collider 

f~i ' (LHC) at CERN Q. Furthermore, in the core of the neutron stars, cold but high-density quark matter with a color 

OhI superconductivity may be realized as has been extensively discussed in recent years \^. 

f^ In relation to QGP, it is of great importance to theoretically understand the phase structure of QCD at finite 

qj [ temperature and density |^. The lattice QCD simulation is one of the most powerful methods for this purpose 
f^ . since it is based on the first principle QCD. Although the finite density simulation is still not available, the phase 
s/ ' transition at finite T in full QCD is extensively studied on the lattice |j]. Another method, which could be potentially 
• ^H , very useful to study the QCD phase structure, is the Schwinger-Dyson (SD) approach |^^. This method is based 
^^ • on the continuum field theory and has been applied for analysing the spontaneous chiral symmetry breaking and 
confinement in the vacuum. Exact SD equations are infinitely coupled integral equations and their solutions yield the 
n-point Euclidean Green's functions. In the practical applications, one needs to truncate the coupled SD equations 
by introducing some ansatze. 

When one studies the spontaneous chiral symmetry breaking at zero temperature and density in the SD approach, 
one starts with the dressed quark propagator in the Euclidean space ||^ , 

S{p) = 



= -i^aA{p')+(TB{p'). (1.1) 

Here, A{p'^) is the wave-function renormalization and B{p^) is the quark mass-function. If one finds a solution 
B{p^) 7^ in the chiral limit, it implies the spontaneous breaking of the chiral symmetry [g]. If the free gluon 
propagator in the Landa u-gau ge and the bare quark-gluon vertex are used in the 1-loop SD equation for the quark, 
A{p^) = 1 holds exactly y,|l^. However, this is not the case for other models of the gluon propagator and/or the 
quark-gluon vertex. The SD equation for the quark is obtained by the extremum condition of the effective potential 
under the variation of the dressed quark propagator S |ll|| . There are two distinct solutions of S in the chiral limit 
corresponding to the Wigner phase and the Nambu-Goldstone phase. The true ground state should be determined 
by calculating the effective potential in those two phases. 

The introduction of temperature or baryon-density to Euclidean QCD reduces the 0(4) symmetry to 0(3). As a 
consequence, the dressed quark propagator in the imaginary-time formalism has a general form [p[ , 
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= -if-paAip'^,^1) -^4t^„C^c(P^W,2J+C^B(P^W^)• (1-2) 

Here a;„ = (2n + l)7rT + i^(n G Z) is the Matsubara-frequency for the quark. T and /i are the temperature and the 
quark chemical potential respectively. j4(p^,a;^) and C{p'^,Lu'^) represent the quark wave-function renormalization, 
and B{p'^,oj'^) is the quark mass-function. For T ^ and/or fi j^ 0, A{p'^,ujl) = C{p'^,uj'^J = 1 does NOT hold even 
if the free gluon propagator in the Landau-gauge and the bare quark-gluon vertex are used in the 1-loop SD equation. 

The purpose of this paper is to obtain a dressed quark propagator at T 7^ by solving the 1-loop SD equation for 
the quark without any ansatz for the functions A, C and B, and to investigate their physical implications. The exact 
numerical solutions of A, C and B have explicit {p'^,uj'^) dependences. Although such functions have been obtained 
in a simple model where the SD equation can be solved algebraically M, the (p^,w^) dependence of A and C for 
realistic models such as QCD-like theory [Q has not been known. Non-trivial (p^jW^) dependence of the functions 
A and C to the nonperturbative phenomena of QCD in the Wigner phase will be also discussed. 

The paper is organized as follows. In Sec. II. A, the SD equation and the CJT effective potential for the quark 
are summarized. Then a model for the gluon propagator, the quark-gluon vertex and the QCD coupling constant is 
introduced at T = 0. In Sec. II. B, they are extended to T ^ 0. In Sec. Ill, basic parameters of the model introduced 
in Sec. II are determined by the condition fj^ = 93 MeV at T = 0. Then the quark wave-function renormalization 
A(p^,w^) and C(p^,Ci;^J and the quark mass-function B{p'^,uo'^) are obtained self-consistently at T ^ by solving 
the SD equation with simultaneous iteration. By calculating the effective potential corresponding to those solutions, 
the true ground state is determined at given T. The critical temperature T^ at which the chiral symmetry is restored 
is also obtained. The effect of the (p^,a;^) dependence of A and C on the chiral phase transition is examined. In 
Sec. IV conclusions and discussions are given. In Appendix, the explicit forms of the integral kernels which appear in 
the SD equation for the quark are given. In this paper, massless (u,d)-quarks are considered, and all calculations are 
done in the imaginary-time formalism |l3]. 

II. SD EQUATION AND EFFECTIVE POTENTIAL 

A. SD eq. for the quark at T = in the rainbow^ approximation and the CJT effective potential 

At T = 0, the SD equation for the quark in the chiral limit is expressed as, 

S-'{p)^S-^^{p) + np)- (2.1) 

Here, Sijj) and S(Q\{p) are the dressed quark propagator in eq. ( |l.l| ) a nd t he free quark propagator l/iji in the chiral 
limit, respectively. S(p^) represents the self-energy of the quark. Eq. ( p.l| ) is derived by the extremum condition for 
the CJT effective potential V[S] under the variation of S{p) |1^; 

V[S] = V^[S] + V2[S], (2.2) 

Vi[S] = j ^tr {in [s-^^{j>)S{p)] - Sl^^{j>)S{p) + l} , (2.3) 

V2[S] = -^J^ {tr [S(p)5(p)]} . (2.4) 

Here "tr" is taken over the Dirac, flavor and color matrices. Vi [S] corresponds to the 1-loop potential with the quark 
1-loop diagram and V2[S] is the 2-loop potential with the one gluon exchange. S(p^) has a general form; 

S(p) = J 0^9Y^,Sik')9K(.P, k)D,Ap - k), (2.5) 

where T'^{p, k) and D^^(p — k) are the dressed quark-gluon vertex function and the dressed gluon propagator respec- 
tively. 

Eq. (2.1) can be decomposed into two coupled integral equations for A{p'^) and B{p^) by taking the trace "Tr" 
over the gamma matrices ||5|,|ljl ; 



A{p') = l-j^Tr\^J:{p)], (2.6) 

B(p2)^iTr[E(p)]. (2.7) 



Since S(p) contains the dressed quark propagator, right-hand sides of eqs. (2.6) and (2.7) are impHcit functions of 
A{k'^) and B{k^) integrated over the momentum k. Once one obtains the solutions Ajjp^) and B{p^), the effective 
potential corresponding to those solutions are given by 

V[S] = -AN^j -0^^ {In (/A2(/) + S2(/)) +pV^(/) - In/ - l} 
= -^ r dp, Hdlplp^ ln{ip^+pl)A^ip^,pl) + B^ip',pl)) 

"" J-oc Jo '^ 

HP^+pl)cTAip^pl)-Hp^+pl)~l }. (2.8) 

To make the connection to T ^ case clea r, th e integration over the 4- momentum are decomposed into \p\ and p, 
instead of p^ = P^ + P, in eg. (2.8) . In eq. ( |2.8| ) the factor 2 for the degenerate (u, d)-quark is included explicitly. 

In order to solve eqs. (2_^) and (2/7) in terms of A{p'^,p1) and B{p'^,p1), one needs to define the structure of the 
dressed gluon propagator D^j/, the dressed quark-gluon vertex funciton F'^ and the coupling constant g^. In this 
paper we adopt the following model: For the dressed quark-gluon vertex function, we use the rainbow approximation 
F°(p, k) = ^7/i, where A" is the color Gell-Mann matrices. For the coupling constant g^, the QCD running coupling 
constant in the one- loop order with an infrared regulator pc is used to describe the momentum dependence iTS] ; 

47r 9 ln[(p2+p2)/AQCDT 

Pc approximately divides the momentum scale into the infrared region and the ultraviolet region. Furthermore, we 
adopt the Higashijima-Miransky approximation |l6yi7|] for the QCD running coupling constant 



g^ Up - kf) =. 9{p^ - k')g'{p^) + e{e p')g\e) = g^ (max[p2, fc^j) (2.IO) 

In this approximation the quark mass- function B{p^) in the ultraviolet region has an asymptotic form consistent with 
that obtained by the operator product expansion and the renormalization group. For the dressed gluon propagator, 
we use the free one 0]^^ (p — k) in the Landau-gauge, 

^1.°.Hp) = (^.. - ^) ^ . (2.11) 

The Landau-gauge is a fixed point of the renormalization group and the gauge parameter is zero to all orders in 
perturbation theory [n8| . 

The integration over the angle between p and k can be performed in (2.5) because the running coupling constant 
has no angle dependence. Then the quark self-energy has no term with odd number of 7^, that is Tr[j^I](p)] = 0. As 
a consequence, eq. (^^) reduces to A{p'^) = 1 in this model as mentioned in Sec. I. The SD equation for the quark 
becomes a single integral equation (|2.7|) for B{p^) after integrating over the angle between p and fc. 



Rr„2. 3Cp f g^umx^l^ k^Bjk') 

^^ ' 16n^ J max[p2,P] -^+^2(^2)' l^-^^i 

where Cp — {N^ — i)/'2Nc is the quadratic Casimir operator for the color SU(A*'c) group. Eq. ( |2.12|) can be expressed 
in term of the variables p and p4 instead of p'^ — p'^ + pi, 

B{p\pl) = j_Jk,j^ d\k\ ^,^^,l^,^^,^^,y E,{\plp,-\k\M). (2.13) 

The explicit form of £'i(|p|,p4; |fc|,fc4) is given in Appendix. Eq. (2.13) has two distinct solutions; one is a trivial 
solution Bw{p^,pi) = representing a free quark, and another is a non-trivial solution Bng{p'^,P4) ¥" 0- The 
subscprit W represents the "Wigner" solution describing a phase in which the chiral symmetry is not broken. The 
subscprit NG represents the "Nambu-Goldstone" (NG) solution which represents a phase in which the chiral symmetry 
is broken spontaneously. To determine the ture gound state realized at T = 0, the effective potential V[S] with the 
Wigner solution and that with the NG solution have to be compared. 



B. SD eq. for the quark at T / 



In this subsection, we extend the formulas in the previous subsection to T ^ 0. For this purpose, we apply the 
imaginary time formalism [nSl and make a replacement 






(2.14) 



where uj^ = (2?^ + l)7rT (n e Z) is the Matsubara-frequency for the fermion. The SD equation for the quark (2.1 
T 7^ is decomposed into three coupled integral equations for A(p^,cj^), C(p^,cj^j) and B{p'^,uj'^) 
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(2.15) 
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+ 2C{k'^,Ujl^)uJjnPjDf^ {P -k,UJ„- LOm) | 



(2.16) 
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(2.17) 



By making the integration over the angle between p and fc, eqs. ( 2.15 ), ( 2.16 ) and ( 2.17 ) reduce to the following 
equations: 



B{p\uI)^2t:T J2 h\^ 



k^B{k,Lu 



2 , ,2 



k'A^k^uD+ulCHk^uD+B^k'^Lol) 



Eii\p\,UJn;\k\,Um); 



(2.18) 



■• m ^ — rv-i \ " 



m— — OO 



{k',col)+ujlC^k^iol) + B^{k\u;l) 

x{A{k^,Ujl^)E2{\p\,UJn;\k\,UJm)+0JmC{k^,0jl,)E3{\p\,UJn;\k\,UJ^)} , (2.19) 

fc2 



^" mf^oo^ fc'A2(fc^O+a;^C2(fc^C^^)+S2(fe^^^) 

X {A(fc^w^)^4(|p|,u;„;|fc|,w,^)+w„C(fc^w^)^5(|p|,w„;|fc|,w,„)} 



(2.20) 



Explicit forms of i?i(|p|, a;„; |fc|, a;,„) {i = 1, 2, 3, 4, 5) are given in Appendix. Eq. (^.18 ) may have two distinct solutions: 
B = Bwip^TUj"^) = and B = BNciP^^^n) ¥" 0- We call the solutions of eqs. (2.1£) and ( 2.20| ) corresponding to Bw 
as {Aw,Cw) and to Bng as {AmgtCng)- The dressed quark propagators in the Wigner phase and the NG phase 
are defined by 



Sw{P,t^n) 
SNG{P,^n) 



i-f-pAw{p^,i^l) +i"fiUJnCw{p'^,(^n) 

1 

~ ij ■pANG{p^,i^l) + i^iUJuCNcip^ ,ojn) + Bng{p'^,ojI) 

= -ij-paAj^oip'^^^l) - *74WnC^C«G(P^'^n) + <^ B jv g ip'^ ^ ^n) 



(2.21) 

(2.22) 



Unlike the case at T = 0, the quark self-energy E(p^, a;,";) can have terms with odd powers of 7^ at T 7^ 0. Therefore, 
A and C cannot be equal to unity even when we take the rainbow approximation with the free giuon propagator in 
the Landau-gauge. Furthermore, A and C have non-trivial dependence on p^ and w^. n 



Once one obtains A,C and i?, the effective potential at T / as a direct generalization of eq. (2.14) can be 
calculated through the formula, 

V[S] = -^T £ fd\p\p' {ln{p'A'ip^col) + o.lC'{p',.l) + B^p^.l)) 

n— — 00 

+p'aA{p\u;l) + u;lac{p\col) - In (p^ + ^2j _ ^ | ^ (2.23) 

To study the chiral phase transition, the difference of the free energy between the NG phase {Bng 7^ 0) and the 
Wigner phase {Bw — 0) should be considered: 



y(r) = v[Sng] ~ v[Sw 

2N, 



^T 5: Jd\p\p^ {in 



n— — OQ 



P^A%aiP^<) + <C%aiP.<) + BlGiP^<) 



p^Al^ip^u^D+ujlC^^ip^cul) 
+P^ [o•A„G(P^Wn) -0-A,^,(p^W^)] +UjI [crc„G(P^W«) -o•cw(p^w^) 
If V{T) > 0(< 0), the chiral symmetry is restored (spontaneously broken). 



}■ 



(2.24) 



III. CHIRAL PHASE TRANSITION AT FINITE TEMPERATURE 
A. Determination of the parameters at T = 

In this subsecti on, w e determine the parameters of the model in Sec. II. A at T = 0. First of all, the momentum 
integrals in eqs. (2^) and ( 2.13| ) are regulated by the ultraviolet cutoff A. As far as A is large enough, physical 
quantities such as /tt and V[S] do not depend on A, while the A-dependence of the chiral condensate {qq)A is 
governed by the renormalization group equation. In our numerical calculation, we adopt A = 5.42 GeV above which 
/tt and V[S] are insensitive to A. 

The chiral condensate {qq)A is known to be insensitive to the infrared regularization parameter pc ||l2|. Therefore 



we takep^/AqcD 



oO.i 



and determine Aqcd to reproduce the pion decay constant f-^ = 93 MeV ]13] . As mentioned 



in Sec. II. A, eq. ( ^.13 ) has two distinct solutions. One is a trivial solution _Bvi/(p^,p|) = corresponding to the free 
quark. For this solution, /tt = and {qq)A — 0. Another solution is B mgJP ^ tPa) 7^ which represents the spontaneous 
chiral symmetry breaking. This solution is obtained by solving eq. ( 2.13| ) numerically by iteration, starting from an 
arbitrary given trial function BtriaKp^iPl) ^ 0. In this case, /^ can be calculated from the formula [[l9[ . 



i^^^ = :;^ ^^4 / d\p\ p^B{p^,pI) I crA(JB + gP^ {cr'ACTB + cja<j'b) \ , 



(3.1) 



with a'g = d(JB{p'^,p'i)/dp'^ and cr^ = daA{p'^,p'i)/dp'^. N^ is a canonical normalization constant for the Bethe- 
Salpeter amplitude in the ladder approximation. 



Nf 



2 
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"■ J-A Jo [ 



tP 



p' 



{aAo"A {cj'^f) +pI {oco'i - {<y'cf) + ^bo'^ {o'^f 



(3.2) 



'^In a simple algebraically soluble model, it was noted before that this non-trivial dependence has significant effect on the bulk 
thermodynamic quantities H. 



with a'c = dac{p^,pl)/dp'^ = a'^ at T = 0. If A = C = I, one obtains N^ = U ^l and Aqcd = 734 MeV is 
determined by the input /^ = 93 MeV in eqs. ( p^ ) or (3^) with A = C = 1 and B^c- Our parameters are consitent 
with those obtained in p^. 



The effective potential for the solution B^a is obtained from eq. (2.8). The difference of the vacuum energy between 
the NG phase and the Wigner phase reads V{0) = —(172.7 MeV)'' < 0. This means that the state corresponding to 
the NG solution is realized at T = 0. At the scale A, the chiral condensate {qq)A is calculated as 

{qq)A = - / 7^:;;^*'^['^(p)] 

o /-A ^A 

= t/ dp4 d\p\p'^aB{p'^,pl). (3.3) 

^ J -A Jo 

With the parameter set determined above, we obtain {qq)A = (—293 MeV)"^. If we change the scale from A to 1 GeV 
by the perturbative renormalization group equation, (qq) igcY — (—217 MeV)^ is obtained, which is consistent with 
the known phenomenological value. 

B. The exact numerical solutions of the SD eq. in the rainbow approximation at T 7^ 

At T ^ 0, we need to replace the dp4 integral by the Matsubara sum, 

/A 7 "max 

^/(P,P4)-T J2 /(P''^"), (3.4) 

'-^ n=-nmax 

with n„ax being the largest integer n satisfying tj„ — (2n + 1)ttT < A. In order to study how the quark wave-function 
renormalization affects the nature of the chiral phase transition, eqs. ( 2.1^ ), ( 2.19| ) and ( p.20D are solved in two 
different cases with the parameter set determined at T = 0. 



(I) : A and C are assumed to be 1 and only eq. ( 2.18| ) for B is solved numerically. The trivial solution i? = of 



eq. (2.18) corresponds to the free quark and the numerical solution Bng 7^ corresponds to the phase of the 
spontaneous chiral symmetry breaking. We define the difference of the free energy between the NG phase and 
the Wigner phase in this case as V(^j-){T) according to eq. (2.24) with Aw = Cw = A^g = Cng = 1- 



(II) : Three coupled eqs. ( ^.18 ), ( 2.19 ) and ( |2.20|) are solved numerically with initial trial functions for A,B and 



C. We call the solutions corresponding to the Wigner phase as Bw{= 0),^vf(7^ 1) and Cwi"^ 1) while those 
corresponding to the NG phase as BMci"^ 0),^jvg(7^ 1) and Cng{¥^ !)■ The difference of the free energy 
between the NG phase and the Wigner phase is defined as V(//)(T) according to eq. ( p.24| ). 

Fig. 1 shows the momentum dependence of the quark mass-function _Bjvg(p^,Wq) in the case (I) with a zeroth 
Matsubara-frequency wq = ttT at T = 10, 150, 200 and 216 MeV. BMciP^ ^^n) i'^ the low momentum region decreases 
with the increase of \n\ at given T. For T > 217 MeV, eq. ( ^.18 ) have only the trivial solution Bw — 0. 



Fig. 2 (a) and (b) show the momentum dependence of the functions Aw and Cw in the case (II) with uoq — nT. 
It is found that Aw and Cw in the low momentum region do not reach unity even at T ~ a few hundred MeV. 
This implies that the quarks and gluons are strongly interacting even in the high T chirally symmetric phase. For 
'-^nin 7^ 0), the difference of Aw and Cw from unity decreases with the increase of \n\ at given T. 

The momentum dependence of the quark mass- function Bno in the case (II) w ith ujq — ttT is shown in Fig. 3. The 



non-trivial solution Bng 7^ is obtained only below T = 155 MeV and eq. (2.18) have only a trivial solution Bw — 
for T > 155 MeV. Bnq for a;„(n ^ 0) in the low momentum region becomes smaller as \n\ increases at given T. 

Fig. 4 (a) and (b) show the momentum dependence of the functions Anq and Cng with ujq = irT. It is found that 
the difference of Ang and Cng from unity in the low momentum region is maximal around T = 150 MeV, and then 
Ang and Cng coincide with Aw and Cw for T > 155 MeV. For tj„(n ^ 0), the difference of Ang and Cng from 
unity decreases with the increase of |n| at given T. 

As mentioned before, in order to determine the critical temperature Tc of the chiral symmetry restoration, o ne ha s 
to calculate the difference of the effective potential between the NG phase and the Wigner phase, V{T) in eq. ( ^.24 ) 



Tc is defined by the temperature at which V{Tc) = 0. Fig. 5 shows the T-dependence of V(^i-){T) and y(//)(T). It is 
found that Tc = 217 MeV in the case (I) and Tc — 155 MeV in the case (II), and the chiral phase transition is of 
the second order in both cases. Due to the crude ansatz yl = C = latT7^0in the case (I), V(/)(T) has a strange 



behavior in which it decreases first and then increases as T increases. Such behavior disappears in V(//) (T) where A 
and C are solved together with B. V(//)(T) increases monotonical ly a s T increases and vanishes at Tc 



The pion decay constant /^ at T 7^ is obtained from eqs. ( |3.l| ) and (^) using eq. (|3.4|). Fig. 6 shows the 
T-dependence of /^ for the NG solution in the case (I) and (II). In the case (I), f-^ has a strange behavior in which 
it increases first and then decreases as T increases because of the crude approximation A — C — I. This behavior 
is correlated with that of V(/)(r) in Fig. 5. On the other hand, /^ in the case (II) decreases slowly first an d then 
changes suddenly near Tc as T increases. The chiral condensate (qq)A at the scale A are calculated by eq.(|3.3|) using 



eq.(3^) at r 7^ 0. In Fig. 7, the T-dependence of {qq)A for the NG solution in the case (I) and (II) is shown. For the 
same reason as that in /tt, the strange behavior of (99) a in the case (I) is found, while {qq)A in the case (II) has an 
almost constant value first and then decreases suddenly near Tc as T increases. 

As mentioned above, the phase transtions in (I) and in (II) are both of the second order, so the critical exponents 
can be defined. Fig. 8 shows the scaling behaviour of {qq)A h^ (I) and in (II) near the critical temperature in each 
case. The critical exponent f3 is defined by 

{qq)A ^ (1 - f ) for T ^ T, - . (3.5) 

(3 is determined by using the linear-log fit numerically as follows 

ln{qq)A = /3 In fl - — j + const. , (3.6) 

where const, is independent of T. We find that (3 ^ 0.51 in (I) and r^ 0.53 in (II) respectively. Both values are 
consistent with (3 in the mean-field theory and agrees with ref . ||l3,E0| . The mean- field critical exponents are a result 
of the structure of the gap equation for the fermion. The inclusion of the soft collective modes will move the mean 
field exponents to the 0(4) values as expected from the renormalization group argument P]. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have obtained the dressed quark propagator as the numerical solutions by solving the SD eq. for 
the quark without any ansatz and studied the effect of the quark wave-function renormalization on the chiral phase 
transition at T ^ 0. The SD eq. for the quark has been solved at one-loop level. We use the free gluon propagator 
in the Landau-gauge, the bare quark-gluon vertex, the one-loop QCD running coupling with the infrared regulator 
and Higashijima-Miransky approximation. The three-momentum and Matsubara-frequency dependent functions for 
the quark wave-funcition renormalization differ from unity especially in the low momentum region. As a result the 
critical temperature decreases substantially compared to that with A = C = 1. 

To examine the effect of the quark wave-function renormalization, we have solved the SD eq. for the quark in two 
cases: 

(I) : the SD eq. for B is solved with an assumption A = C = 1. 

(II) : coupled SD eq. for A, B and C are solved numerically. 

(I) is an approximate solution, and (II) is the exact solutions of the SD eq. for the quark. The Wigner solution in the 
case (I) corresponds to the free quark. From the difference of the free energy between the NG phase and the Wigner 
phase in the case (I), the critical temperature of the chiral phase transition turns out to be Tc = 217 MeV. Also 
the difference of the free energy between the NG phase and the Wigner phase in the case (II) shows Tc = 155 MeV. 
Exact numerical solution with A^ 1 ^ C in the case (II) not only reduces the critical temperature but also removes 
the pathological behavior of /^r and {qq) in the case (I) . We found that after the chiral symmetry restoration, quarks 
are not free even at temperature of a few hundred MeV in the case (II). In both cases (I) and (II), the chiral phase 
transition is of the second order. The critical exponent (3 is extracted from {qq) near Tc and we obtained (3 ^ 0.51 in 
(I) and /?-- 0.53 in (II). They are consistent with the result of ||,||J20|. 

Finally, several comments are in order. The gluon propagator used in this paper is the free one and cannot explain 
the quark and gluon confinement. To study the effect of the quark confinement on the chiral phase transition, it would 
be necessary to use the phenomenological gluon propagator such as the dual Ginzburg-Landau model which explains 
the spontaneous chiral symmetry breaking and the quark confinement pl^,E2| . The screening of the gluon should be 
also included. The electric sector of the gluon propagator is infiuenced by the Debye screening and the magnetic sector 
by the non-perturbative magnetic mass in the static case or the dynamical screening in the non static case. The study 



of the chiral phase transition at finite density and temperature is also an important future problem in our approach 
to locate the triciritical point in the QCD phase diagram [^ . Our preliminary study shows that the location of the 
tricritical point in our model is {Ttr,n^^) = (210 MeV, 129 MeV) for the case (I) and (142 MeV, 246 MeV) for the 
case (II), where /x^(= 3/Xtr) is the baryon chemical potential. 
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APPENDIX: 



In this a ppend ix, we show explicit forms of five kernels i?i(|p|,cj„; \k\,ujm) {i — 1,2,3,4,5) in eqs. (2.7), (2.18) 



(2.19) and (2.20) 



4 1 

-Bi(|p|,Lj„; \k\,u;,n) = -- 



97r|p| |fc| In [(max[p2 +ujlk^+ cul] + pi) /Aqcd'] 

xin ;i^|-|;i;; +;—-;; (ad 

{\P\ + \k\)^ + (uJn ~ C^m)2 



Q 1 

E2{\p\,ujn; \k\,ujm) = -;;;;- 



277r In [(max[p2 + ojlk^ + col] + p^) /Aqcd'] 

j g P^ + fc^ + 3{uJn - UJynf , (IpI - \k\f + (CJ„ -i^mf 

""^ 4|p||fc| ■ "^ (IpI + |fc|)2 + (c.„ - .;„)2 

(p2 - fc^ + (CJ„ - ^ra?) (-p2 + fc^ + (CJ„ - UJ^^f) 

4|p||fe| 

1 1 



-B3(|p|,Lj„; |fc|,cj„j) = - 



(IpI - |fc|)2 + (CJ„ - i^ra? (IpI + \k\f + (CJ„ - U.^Y 
4(^„ - LJ„t) 1 

277r|p||fc| ■ In [(max[p2 + ,,2^ fc^ + ,,2j +^2) /Aqcd'] 



^ J 1 (|p|-|fc|)^ + (^»-^m)^ I / „2 , . 2 , / ^2^ 



1 1 



i?4(lp|,W„; |fc|,Wm) = 



(IpI - |fe|)2 + (tj„ - tj„,)^ (|p| + |fc|)^ + (^« - ^"0^ 

4(u;n - u;m) 1 

27^|p||fe| ■ In [(max[p2 + ^2^ fe^ + u^ + p2) /Aqcd'] 

^ /, (|P|-|fc|)^ + (^n-^m)^ , / 2 .2 , / ^2^ 

n'"(iPi+ifci)^+(-.— ™)^''^^ "' +K--™)) 
1 1 



(IpI - |fe|)2 + (tj„ - 1^„,)2 (IpI + |fc|)2 _^ (^^^ _ ^^^j2 

4 1 

i?5(|p|,-„; |fc|,-„0 = ^^^^^ • l,[(^^,[p2+^2^fc2_^^2^j^p2)/AQc^2] 
r (|p|-|fc|)2 + (^„-^„)2 

y (IpI + |fc|)2 + (cj„ - wm)2 

1 1 



(IpI - |fc|)2 + (a)„ - 1^„,)2 (IpI + |fc|)2 ^ (^^^ _ ^^^j2 
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FIG. 1. The momentum dependence of the quark mass-function Bng{p ,(^o) in ^hs case (I) at T =10,150,200 and 216 MeV 
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FIG. 3. The momentum dependence of the quark mass-function Bng{p ,(^a) in the case (II) at T =10,100,150 and 154 MeV 
with the lowest Matsubara-frequency loq = t^T. 
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